Non-perturbative renormalization-group approach to lattice models 
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The non-perturbative renormalization-group approach is extended to lattice models, considering 
as an example a (f) 4 theory defined on a d-dimensional hypercubic lattice. Within a simple approxi- 
mation for the effective action, we solve the flow equations and obtain the renormalized dispersion 
e(q) over the whole Brillouin zone of the reciprocal lattice. In the long-distance limit, where the lat- 
tice does not matter any more, we reproduce the usual flow equations of the continuum model. We 
show how the numerical solution of the flow equations can be simplified by expanding the dispersion 
in a finite number of circular harmonics. 
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I. INTRODUCTION 

The non-perturbative renormalization group (NPRG) 
has proven to be a powerful tool in the study of systems 
with a large number of interacting degrees of freedom 
It has been successfully applied in many areas of physics, 
condensed matter, nuclear and particle physics, etc. (For 
a review and a pedagogical introduction, see Refs. and 
&) 

The NPRG is based on an exact flow equation sat- 
isfied by the effective action T[M] (i.e. the generating 
functional of one-particle irreducible vertices). In the 
most common approximation, one expands the effective 
action in powers of derivatives of the field M . While 
this approach is often the easiest way to solve the NPRG 
equations, it yields only the small momentum behavior 
of the vertices £ As shown recently, it is possible to solve 
the NPRG equations beyond the derivative expansion 
and compute the whole momentum dependence of the 
vertices, 7 - 8 - 9 - 10 ! 11 ! 12 - 13 - 14 

A proper description of the momentum dependence of 
the vertices opens up the possibility to study lattice mod- 
els. While the derivative expansion is always appropriate 
to study the long-distance physics and predicts universal 
quantities such as the critical exponents, it often fails to 
predict non-universal quantities (e.g. the critical temper- 
ature of a phase transition) which depend on the short- 
distance physics. Moreover, in some cases the lattice is 
the very reason for the occurrence of a phase transition.— 
(For approximate treatments of lattice models within the 
derivative expansion, see Refs. fl^fl7lfl8lfl9l .) 

In this paper, we show how the usual NPRG approach 
should be modified in the case of lattice models. As 
an example, we consider a 4> A theory defined on a d- 
dimensional hypercubic lattice. We solve the flow equa- 
tions within a simple approximation for the effective ac- 
tion. From the 2-leg vertex, we deduce the renormalized 



"dispersion" (i.e. the kinetic energy in particle language) 
over the whole Brillouin zone. In order to simplify the nu- 
merical solution of the flow equations, we introduce two 
approximations: i) the LPA' (where LPA stands for Local 
Potential Approximation) - a natural generalization to 
the lattice case of the LPA' used in continuum models^ 
-, which neglects the renormalization of the dispersion 
except at very small momenta; ii) the H-LPA' where the 
LPA' is supplemented by a circular harmonic expansion 
of the 2-leg vertex, thus allowing one to take into account 
the renormalization of the dispersion in a (numerically) 
efficient way. By comparing with the exact solution of 
the flow equations, these two approximations are found 
to be remarkably accurate. 



II. NPRG FOR LATTICE MODELS 



We consider a 
percubic lattice, 



General method 



theory defined on d-dimensional hy- 



S[<f>] = £ {^OHV) + v]cf> r + -^j 



(1) 



where {r} denotes the sites of the lattice. For simplic- 
ity, we consider a one-component real field r ; the ex- 
tension of our approach to a iV-component field model 
with O(N) symmetry is straightforward. The bare dis- 
persion eo(q) is chosen such that eo(q = 0) = and 
liniq^o eo(c[) = eoQ 2 - (This is always possible by a re- 
definition of the parameter v in (Q]).) For a system with 
nearest-neighbor interactions only (the case we shall con- 
sider for the numerical solution of the flow equations) , 



£ o(q) 



d 

-2e £[cos(&,) - 1], 



(2) 
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where the maximum of eo(q) is given by 4cfeo- The lattice 
spacing is taken as the unit length. 

To implement the RG procedure, we add to the action 
the regulator term 



\ ]T r i?HV)0 r = \ 0- qJ R(q)0 q , (3) 



where the function i?(q) = i?/c(q) depends on the energy 
q = N~ x / 2 ^ r e~ lq ' r </> r is the Fourier transformed 
field with N the number of lattice sites. The sum over q 
in ([3]) is restricted to the first Brillouin zone ] — 7r,7r] d of 
the reciprocal lattice. In the thermodynamic limit (N — > 
oo), 



i v - r ^ r ^ 

AT 2^ ^ _/ 2 7T " ' 2^ 



(4) 



The NPRG approach is based on the effective action 
r[M] defined as the Legendre transform of the free en- 
ergy — lnZ[/i] from which ASr[M] is subtracted. Here 
h denotes an external field that couples linearly to the 
4> field and M T = {4> r )h is the expectation value of 4> r . 
T[M] satisfies the exact flow equation^. 

d k F[M] = ±Tr[d k R(rW [ M ] +RY 1 } (5) 

as the energy scale tk is varied. r^ 2 '[M] is the second- 
derivative of T[M) with respect to M. To keep the no- 
tation simple, we do not explicitly indicate the fc depen- 
dence of T, r( 2 ) and R. 

An important difference with continuum models comes 
from the regulator function R. The latter is chosen to be 
of the form 



R(q) = Zeo(q)r(e (q)), ?o(q) = 



eo(q) 



(6) 



where the fc-dependent variable Z will be specified below. 
A typical choice for the function r is r(x) — (e x — 
The regulator function R gives a mass of order efe to 
low-energy fluctuation modes (eo(q) J; e fc) but leaves the 
high-energy modes (eo(q) ^ essentially unaffected. 
When tk —>■ oo - or, in practice, tk larger than any typical 
energy scale - all fluctuations are frozen and mean-field 
theory becomes exact: T[M] = S[M], As long as tk 3> eo 
(i.e. efc ^> eo(q)), fluctuations are local (on-site). They 
become non-local when ~ eo- Although i?(q) = i?fc(q) 
is a function of ek, it is convenient to write €k = eok in 
terms of a momentum scale fc and consider all quantities 
of interest to be function of k rather than e^. The regime 
where fluctuations are local in space then corresponds to 
length scales fc" 1 <C 1, i.e. much smaller than the lat- 
tice spacing, whereas the condition fc -1 ^ 1 implies that 
fluctuations can propagate through the lattice. Since the 
function r(x) typically vanishes exponentially for x >• 1, 



!imi* ( q) = Wr(^ 



(7) 



and we reproduce the regulator function that is used in 
continuum models (eo(q) = £oq 2 )-— ^ Eq. J7]) expresses 
the fact that when £k <C to (i.e. k <c 1) the lattice does 
not matter any more and only the small-q limit eoq 2 of 
the dispersion plays a role. To distinguish between these 
two regimes, characterized by the presence or absence of 
strong lattice effects, it is convenient to introduce the 
crossover momentum scale k x ~ 1 (e^ ~ eo)- 

In the following we write k = k(t) = Ae* (t < 0) where 
A is such that when = eoA 2 , all fluctuations are effec- 
tively frozen and the mean-field theory a good approxi- 
mation. In practice, one should verify that the solution 
of the flow equations remains essentially unchanged when 
A is increased above the chosen value. 



B. Flow equations 

In this section, we derive the flow equations in the case 
where the effective action is approximated by the simple 
form 



r[M] 



r 



)M T + U( Pr ) 



(8) 



where e(q) denotes the (fc-dependent) dispersion. For 
symmetry reasons, the potential U can only be a function 
of p r = M 2 /2. In the ordered phase - the only one that 
will be of interest to us-, we approximate it by 



U(Pr) = ^(Pr - Po) 2 , 



(9) 



where po is the fc-dependent minimum of the potential. 
The wave-function renormalization factor Z is defined by 



lim 6 q 
q^o e (qj 



lim 



e(qj 



(10) 



The dependence of R on Z [Eq. ((6])] is a necessary con- 
dition for the effective action T[M] to reach a fixed point 
when the system becomes critical (po{k —> 0) = + ). We 
are therefore left with three unknown parameters (po, A 
and Z) and a function e(q) to be determined as a func- 
tion of k{t). The initial values at t = are po = — 3v/u 
(for v < 0), A = u/3 and e(q) = eo(q), where u and v are 
introduced in (Q]). 

Given our choice of the regulator function [Eqs. (07])], 
we expect the flow equations to reproduce those of the 
continuum model when k <C 1. This suggests to intro- 
duce the dimensionless variables 



M r = (Ze k ) 1/2 k- d / 2 M r , 
A = Z- 2 e^ 2 k d X, 



p r = Ze k k p r , 
e(q) = (ZekJ-Mq). 



(11) 
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The effective action then takes the form 

r[M] = fc d ]T 



d d r 



i(v f Af(r)) 2 + ^ (/ 5( r )_p ) 2 



(12) 



where the last result is obtained in the continuum limit 
using limfe-,0 — q 2 /fc 2 and introducing the dimen- 
sionless continuous variable f = kr. 

The flow equations are deduced from and iJH). In 
dimensionless form, one finds (Appendix [A| 



dtPo 
d t ~X 

V 



-(d-2 + rfipo + 3Lf(k, 2\p , m e), 
(d - 4 + 2r;)A + 9A 2 L^(fc, 2Ap , »7! e), 

36A 2 p Af 2 d (fc,2Ap ;e) 
l-36A 2 poMf(fc,2Ap ;e)' 



(13) 



where the threshold functions L d and Mf 2 are defined 
in Appendix [HI and 



-9A 2 p fc" rf / { [(r?r + 2e r')e g 2 ] g(q' + q) 
+ff(q') [(r?r + 2e r')e 5 2 ] q , +q 
-2[(r ? r + 2e r / )?o5 2 ] q ,5(q / )}: 



(14) 



where r = r(efo(q)) and r' = r'(eo(q))- We have intro- 
duced the running anomalous dimension 



i] = -d t In Z 
and the dimensionless propagator 

s(q) ' 



e(q)+e (q)r-(e (q)) + 2Apo 



(15) 



(16) 



A direct consequence of the property |(7j) is that the flow 
equations lfl3|) become identical to those of the contin- 
uum model in the limit k <C 1, 



d t po = -(d-2 + r))p + 6v d lf{2\p ,n), 
d t \ = {d-± + 2 n )\+l%v d \Hi{2\p Q ,rj), 

rj = 72^\ 2 p m d (2~\p , V ), 
d 



(17) 



where i„ and m d are the usual threshold functions (Ap- 
pendix [Bj, whereas 1|14|1 becomes identical to the self- 
energy equation derived in Ref. 0. The coefficient Vd is 
defined in Appendix iBl 



<lv Qv'), the numerical solution of the flow equations 
will be quite demanding in particular for d > 3. In this 
section, we discuss two approximations which make the 
numerical solution much easier. Their reliability will be 
discussed in Sec. IIIII 



1. LPA ' 

The numerical solution of the flow equations lfT3|) and 
(fl4|l gives a complete description of the behavior of the 
4 theory |(T]) in the approximation where the effective 
action T[M] is given by ([8]). It yields not only the crit- 
ical exponents but also the propagator G(q) = (4>q4>-q) 
over the entire Brillouin zone. We know from previous 
works 7 *^^ that the anomalous dimension 77 can be ob- 
tained either from Z [Eq. lfT5|l ] or from the momentum 
dependence of the propagator G(q) ~ l/|q| 2_,? in the 
limit q — > 0. In order to simplify the numerical solution, 
we can give up the exact solution of the flow equation 
(fT4|l satisfied by e(q), in particular near q = 0, since we 
still have the possibility to extract r\ from the scale de- 
pendence of the wave-function renormalization factor Z 
[Eq. mi 

In the simplest approximation, one sets e(q) = eo(q), 
i.e. e(q) = Ze (q). For e% <C eo, this yields e(q) = 
Zeoq 2 , which is nothing but the LPA' previously intro- 
duced for continuum models. The LPA' amounts to writ- 
ing the effective action Y[M] in terms of a local potential 
U (with a possible expansion in field truncated to a finite 
order) and a gradient term Z(VM) 2 whose amplitude is 
given by Z. As it is based on a gradient expansion, the 
LPA' is valid at small momentum |q| < k and gives only 
the long-distance behavior of the propagator. It is made 
possible by the fact that only propagators with momenta 
|q| ^ k enter the flow equations. In the lattice case we are 
considering here, the approximation e(q) = eo(q) can be 
seen as the natural extension of the LPA' used in contin- 
uum models. In the limit ek *C to where the lattice plays 
no role any more, its accuracy will be identical to that of 
the LPA' in continuum models. When 3> eo, fluctua- 
tions are local in space so that the dispersion should not 
be significantly renormalized: Z ~ 1 and e(q) ~ eo(q)- 
Whether or not the LPA' is also reliable when eu ~ eo 
is more difficult to assess without comparing to a more 
complete solution of the flow equations. This will be 
done in Sec. IIIII The flow equations in the LPA' reduce 
to (jT3j) where the threshold functions L d and Mf 2 should 
be computed with the replacement e(q) — > ?o(q). 



C. Approximations 



2. Circular harmonic expansion 



Even if one takes advantage of the symmetry properties 
of the dispersion e(q) (e.g. invariance under q v — * —q v or 



It is possible to go beyond the LPA' and obtain the 
renormalized dispersion over the entire Brillouin zone by 



4 



expanding e(q) in circular harmonics, 

L 

e(q) = ^2 e nm i[cos(nq x ) cos(mq y ) cos(lq z ) - 1], (18) 

n,m,/— 



and retaining only a subset of harmonics - defined in lfT8j) 
by the integer L. In this section, we take d — 3. For a 
system with nearest-neighbor interactions only [Eq. 
the initial conditions for the coefficients e nm i are given 
by 



1 4=0 



+S n ,o8m,i8i,o + S n fiS m ,oSi,i). (19) 



The flow equations for the dimensionless amplitudes 
Inmi = t nm i/{Zek) read (Appendix [A} 



dt^nmi — (v — 2)e nm ; — 18c lml /A 2 p 

x ^(*.2Apo,»j;e)^](2AA ) ;c), (20) 

where the functions H^ 1 ' and il d ( 2 ) are given in Ap- 
pendix El and c nm i = (2 - 5 nt o)(2 - S mfi )(2 - <5/ j0 ). 

When 6k ^> eo , the fluctuations are local and the renor- 
malization of the dispersion e(q) should be negligible. 



When efc 



< 



non-local fluctuations renormalize the 



dispersion and induce harmonics which are not present 
in the bare dispersion eo(q). We expect the harmon- 
ics cos(nq v ) to be generated when n ~ k^ 1 . On the 
other hand, a finite value of e(q) acts as a mass term 
in the propagator lfT6|) . so that we expect the flow of 
e(q) to stop when e(q) ~ 1, i.e. when k ~ |q| for 
k, |q| <g; 1. Thus the highest harmonics cos(nq u ) is of 
order n ~ min(fc _1 , |q| _1 ). We deduce that the harmon- 
ics expansion (fLS)) is valid only if L ^> min(fc , |q| _1 ). 

In the limit k — > 0, any finite truncation (L < oo) 
cannot properly describe the renormalized dispersion for 
|q| -C L^ 1 and therefore the q — * limit of the prop- 
agator. In particular, it will always give a dispersion 
that behaves as q 2 for q — > and a critical propaga- 
tor G(q) ~ 1/q 2 with vanishing anomalous dimension. 
As in the LPA', one should then extract the anomalous 
dimension r\ from Z. 

Finally, we note that we can combine the harmonic ex- 
pansion with the LPA' by replacing e(q) by eo(q) in the 
threshold functions L*, Mf 2 and iJ d ( 1 - 2 ). This approxi- 
mation will be referred to as the H-LPA'. 
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FIG. 1: (Color online) p vs -t in the LPA' ((red) solid 
line) near criticality for X(t = 0) = 100 (d = 3). The black 
squares are obtained from the full solution of the flow equa- 
tions l|13|14p . The ( green) dashed line is obtained from the 
flow equations of the continuum model with the same bound- 
ary conditions at t = —16. 
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FIG. 2: (Color online) Same as Fig.[T]but for A. 
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III. RESULTS AND DISCUSSION 



FIG. 3: (Color online) Same as Fig. [JJbut for rj. 



In this section, we consider a three-dimensional system 
with nearest-neighbor interactions [Eq. J21)] and take eo = 
1. 

Let us first discuss the case of a strong initial value 
X(t — 0) = 100 of the interaction. The flow of po, A 
and f] is shown in figures [U [2] and [31 The initial value 
po(t = 0) ~ 0.1358 is adjusted (to a precision of ~ 10 -16 ) 



so that the system is (nearly) critical, as shown by the 
plateaus observed for \t\ > 12 in Figs. QJ3J 21 We find a 
remarkable agreement between the LPA' and the full so- 
lution of the flow equations (|13jl4[> . which shows that 
the LPA' is a very good approximation for any value of 
efc. Due to the numerical cost, we have only attempted 
to solve l|13|14p for a limited range of t when no approx- 
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imation is made. Although this seems to give a slightly 
better estimate of the anomalous dimension (the exact 
value is close to 0.036), the improvement over the deriva- 
tive expansion is very weak. The main limitation of our 
approach wrt a more accurate calculation of 77 comes 
from the neglecting of the p dependence of the renormal- 
ized dispersion e(q) (in particular of the wave-function 
renormalization - factor Z defined in lfT0|)) and the sim- 
ple truncation of the potential U(p) introduced in ([9]). 

Figures [U [2] and [3] also show the flow obtained for the 
continuum model (eo(q) = q 2 ) with the same boundary 
conditions at t = —16. In practice, we take the final 
values at t = —16 and solve the flow equations backwards 
in "time" (i.e. from t = — 16 to 0) with the replacement 
e o(q) - * q 2 in the threshold functions L d and Mf 2 , i.e. 
L d n -> 2vdi and M(t] + M$ -> m d (see Appendix EJ . 
We can clearly distinguish between a short-distance (or 
high-energy) regime where the lattice effects are strong 
and a long-distance regime (k < k x ~ 1 or \t\ > 3.2) 
where the lattice effects rapidly disappear and become 
undetectable in the limit k -C k x iH 

Given the success of the LPA' for calculating p Q , A 
and 77, it makes sense to consider the H-LPA' to com- 
pute the renormalized dispersion. Figures [4] and [5] show 
the dispersion e(q) as well as the first harmonics e nm i 
(n,m,l < 2). The agreement between the H-LPA' and 
the full solution of the flow equations (|13|14p is again 
remarkable. The amplitudes e nm i decrease very rapidly 
with n,m,l, and only a small number of harmonics is 
required for an accurate description of e(q) except - as 
discussed in Sec. Ill C 21 near q = 0. Note that po and 
A are sensitive to local fluctuations, while e(q) (or e nm i) 
is not: as expected d t e n mi — when \t\ < 2. When 
k r~j k x , the amplitude eioo = eoio = £001 of the first 
harmonics varies with k and higher-order harmonics are 
progressively generated (Figs. [5] and EJ. 

In Fig. [7l we show the flow of po, A and 77 for a weaker 
value of the bare coupling, X(t = 0) = 0.1. To understand 
the dependence of the results on X(t = 0), one should 
introduce the Ginzburg scale k c , which is proportional 
to X(t = 0) in three dimensions . 8 ' 13 ' 14 In continuum 
models, the Ginzburg scale separates the infrared region 
£ "C |q| <C k c (with £ the correlation length) character- 
ized by the scaling form r^ 2 )(q) ~ |q| 2-1 ' of the inverse 
propagator, where rj* = lim^-oo rj is the anomalous di- 
mension, from a perturbative regime k c -C |q| <C Ao, 
The microscopic cutoff Ao should be much larger than 
k c for the perturbative regime to be observable. In the 
regime k c <C k -C Ao, the running anomalous dimension 
is nearly zero and the flow is dominated by the Gaus- 
sian fixed point. Figures CEJ [H and [3 (A(t = 0) = 100) 
correspond to the case where k c >• k x . As soon as the 
lattice scale k x is reached (\t\ ~ 3.2), the flow rapidly 
crosses over to the critical regime characterized by the 
fixed point values p~Q, X* and the anomalous dimension 
77*. Note that the finite value of r)(t — 0) obtained from 
the continuum model (green dashed line in Fig. [3]) is ex- 
plained by the fact that even at t = the system is 
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FIG. 4: (Color online) Renormalized dispersion e(q) obtained 
in the H-LPA' ((red) solid line) and from the full solution of 
|113I14[| (black squares) along the diagonal (0, 0, 0) — > (it, tt, n) 
of the three-dimensional Brillouin zone. The bare dispersion 
£o(q) is shown by the (green) dashed line. 




FIG. 5: (Color online) Harmonic amplitudes e nm i obtained 
in the H-LPA' (lines) and from the full solution of I|13ll4p 
(symbols). At t = 0, only the harmonic eioo = eoio = eooi = 
— 2 is nonzero. 




k 

FIG. 6: (Color online) Harmonic amplitudes e„oo (lines), e„ n o 
(squares) and e nnn (circles) vs k for n = 1,2,3 and 4 (from 
right to left). All amplitudes are normalized to their values 
at k — * 0. As in Fig. \E\ eioo is shifted by 2. 
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FIG. 7: Same as Figs. [TJ [2] and [3] but for the initial condition 
A(t = 0) = 0.1. 



not in the perturbative regime (i.e. k{t = 0) < k c ) for 
\{t = 0) = 100. On the other hand, for X(t = 0) = 0.1 
(Fig. [7l), one has k c -^k x . Between the local fluctuation 
regime and the critical regime, one can clearly observe 
an intermediate regime k c -C k <C k x where the run- 
ning anomalous dimension is nearly zero and the run- 
ning coupling constants po and ^significantly differ from 
their fixed point values p^ and A*. Note that since the 
renormalization of A and r\ is nearly zero in the local fluc- 
tuation regime (k < k x ), the flow of these quantities is 
well approximated by the continuum model equations for 
all values of k (hence the superposition of the (red) solid 
and (green) dashed lines in the corresponding figures). 



IV. CONCLUSION AND PERSPECTIVES 

We have shown how the presence of a lattice can be 
taken into account in the NPRG. Our approach allows 



one to compute both critical exponents and non-universal 
quantities such as the critical temperature or the renor- 
malized dispersion. We have proposed two approxima- 
tions which considerably reduce the numerical difficulty 
of solving the flow equations. While the LPA' is suf- 
ficient to obtain the small-momentum behavior of the 
propagator, the H-LPA' - which is based on a circular 
harmonic expansion of the dispersion - yields the renor- 
malized dispersion over the entire Brillouin zone except 
at very small momenta (as in the LPA', the small momen- 
tum behavior of the propagator is deduced form the scale 
dependence of the wave-function renormalization factor 
Z). It is straightforward to generalize our approach to 
more complicated truncations of the effective action than 
the one considered in this paper [Eq. ([8])] as well as to 
quantum- mechanical many-body systems. 
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APPENDIX A: FLOW EQUATIONS 

The flow equations for the potential U (p) and the dis- 
persion e(q) can be deduced from the flow equation of 
the effective action T[M] and the vertex in a uni- 
form field M r = M = const. Since e(q = 0) = 0, one 
has 



T[M] 

r (2 W) 



M r =M 



M r =M 



NU{p), 

<W,o [<q) + U' + 2pU"] , 



(Al) 



where U' = d p U and U" = d 2 p U . From {5j, one then 
finds 



dtU(p) = ~ j G?(q)fl(q), 



(A2) 



where 



G(q) 



r( 2 >(q,-q)+i?(q) 



e(q) + i?(q) +U' + 2pU" 



(A3) 



is the propagator in a uniform field. Using i?(q) = 
—Zeo(q)[rjr + 2eo(q)r'] and r\ defined by (jT5j) , one ob- 
tains 



dtU = -\j eo(q)[r?r + 2£ (q)r']<7(q) 
= k d L d {k,U' +2pU",r);e), 



(A4) 
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where we have introduced the dimensionless propagator 

1 



g(q) = Ze k G(q) = - 



e(q) + e (q)r + U' + 2pU" 



(A5) 



The lattice threshold function Lq is defined in appendix 
IB! The dimensionless potential U(p) — k~ d U(p) satisfies 

d t U = -dU+{d-2+T))pU'+L$(k, U'+2pU", 77; e). (A6) 

With the truncation ([9]), we obtain the flow equation of 
po from the condition U'(p~o) = 0, i.e. 

d t U'(p ) = d t U' +U"(p )d t po = 0, (A7) 

Po 

which gives the first of equations lfT3|) . Similarly, the flow 
equation of A is derived from A = U"(po), i.e. 



d t X = dtU" 



-U^(po)d t p 



(A8) 



(with =0), which gives the second of equations (fT3| . 
(The flow equations of U' and U" are deduced from (|A6[) 
using (|B2|).) 

The flow equation ([5]) of the effective action implies 
that the vertex r( 2 )(q, q') = (Jq+q'^r^^q) satisfies 

d t T^(q) = ^G(q')r( 4 »(q,-q I q' ! -q') 
q' 

- \dt J2 G(q')G(q' + q)rW (q, q', q' - q) 
q' 

xr (3 »Kq' + q-q') (A9) 

in a uniform field M r — M. d t = (dR/dt)dn acts only 
on the regulator function R. The vertices in (| A9[) are 
obtained by differentiating ([Sj) and setting M r = M = 



r (3) (qi:q2,q3) = £ qi +q 2+q3 ,o 



3Xy/2p~ 



N 
3A 



r (4) (qi,q2,q3,q4) = s^+^+^+^o—. (Aio) 

Since e(q) = r( 2 )(q) - T< 2 )(q = 0), one deduces 



d t e(q) = -9X 2 p d t / [G(q')G(q' + q) - G(q') 2 ] . 
Jq' 

(AH) 

In dimensionless form, we obtain 
9 t e(q) = (r, - 2)e(q) 

- 9~X 2 p k- d d t [ [ 5 (q') 5 (q' + q) - g(q') 2 } , (A12) 
Jq' 

where the propagator g(q) is evaluated at p = po. Using 
d 



dt = R 



dR 



-{ryr + 2e r ) — , 



(A13) 



we obtain (jT4j) . 

The equation for the anomalous dimension is obtained 
by expanding (jAl2|) to order q 2 using ?o(q),e(q) — > 
q 2 /fc 2 for q — * 0, 

r) = -^p k 2 - d d t J [d qi g{ci)f ■ (A14) 

With (|A13|) . we deduce 

t? = -9X 2 p k 2 - d [ [d qai g(q)} 

Jq 

x d q ^ [{ W + 2e~ r')e~o(q).9(q) 2 ] , (A15) 

which gives the last of equations (fTB")) when expressed in 
terms of the threshold functions Mf and M$ defined in 
appendix IbI 

If one expands the dispersion in circular harmonics as 
in lfl8l) . one has 



d t e 



t^-nml — Cnmi 



1 / cos(nq x ) cos(mq y ) cos(lq z )d t e(q) 



-18c nm iZe k X 2 p k d g(r nm i)d t g(r nml ), 



(A16) 



where c nm i is defined in Sec. lIIC"2l and r nm ; denotes the 
position of the lattice site with coordinates (n, m, I). g(r) 
is given by the Fourier transform of the dimensionless 
propagator l|16p . Using (| A13|) to compute dtg(r nm i), one 
eventually obtains l|20p . 



APPENDIX B: LATTICE THRESHOLD 
FUNCTIONS 

1. Definition 

The lattice threshold functions are defined by 

k~ d 



L i(k,w,r);e) = -(n + 6 n0 )- 



x / e~ (q)[r?r + 2eo(q)r'] 5 (q)"+ 1 , 

Jq 

M d (k,w;e) = f [d qm g(q)]d qm [e~o(q)rg(q) 2 ] 

M d (k,w;e) = J[d q ^g{q)]d qsc [? Q {q)r' g[qf] , 

H^](k,w,rr,e) = k~ d [ e^ r -' [ W + 2e~ (q)r'] 

Jq 

x£ (q)g(q) 2 , 

H d n ^(w,i) = / e * qr "™\g(q), (Bl) 

Jq 

where g(q) = l/(e(q) + £o(q)r + w). The functions L d 
satisfy 

d w L*(k,w,r);e) = -(n + 6 ni0 )L* +1 (k,w,T];e). (B2) 



k — > limit where 



In the limit k — ► 0, the functions r and r' ensures that 
the integrals determining L d and Mf 2 are dominated by 
| qj < k -C 1. One can then use ?o(q) — 2/ and e(q) ~ y 
(y = q 2 /fc 2 ) and replace the integral over the Brillouin 
zone by 



dqi 
2ir 



dq q 
2n 



2v d k d / dyy d ' 2 - 
Jo 



(B3) 



where v^ 1 = 2 d+1 Ti d ' 2 T{d/2). This gives 

lim fe ^o L d (k, w, rj; e) = 2v d l d (w, rf), 
lim^o [M d {k, w- e)r) + M d (k, w; e)] = -^m d (w, rj), 



(B4) 



m d (w 1 rj) 



~(n + S n>0 ) dyy d/2 (T 1 r + 2>ir ! )a 



dyy d ^(l + r + yr')g 4 

x{[rjr + (rj + 4)yr' + 2y 2 r"} 
—2yg(l + r + yr )(rjr + 2yr')} 



(B5) 



(g = l/(y(l+r(y))+w)) are the usual threshold functions 
of the continuum model. Note that in the k — ► limit, 
the threshold functions L d and M d 2 become independent 
of k and e(q). 
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